Dynamical breakdown of the Ising spin-glass order under a magnetic field 
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The dynamical magnetic properties of an Ising spin glass Feo.55Mno.45Ti03 are studied under 
various magnetic fields. Having determined the temperature and static field dependent relaxation 
time t{T;H) from ac magnetization measurements under a dc bias field by a general method, we 
first demonstrate that these data provide evidence for a spin-glass (SG) phase transition only in 
zero field. We next argue that the data t{T;H) of finite H can be well interpreted by the droplet 
theory which predicts the absence of a SG phase transition in finite fields. 
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PACS numbers: 75.10.Nr,75.40.Gb,75.50.Lk 

One of the recent interests in the spin-glass (SG) study 
is the nature of Heisenberg spin glasses under a mag- 
netic field P, There exists, however, a still unset- 
tled problem even on the conceptually much simpler Ising 
spin glasses. The mean-field theory predicts that the SG 
phase remains up to the de Almeida-Thouless (AT) line 
in the {H,T) plane The droplet theory, based on 
the short-range Edwards- Anderson (EA) model, instead 
predicts that any applied magnetic field will break up 
the equilibrium SG long-range order 0, 0|- Experimen- 
tal investigations on Ising spin glasses give evidence both 
for the existence of an AT-line Q and against critical 
dynamics under a field 0, while some recent simulation 
and experimental results indicate that the SG phase does 
not exist under a magnetic field [sll^ HolllH . 

In the present paper, we address the problem through 
ac susceptibility measurements on the Ising spin glass 
Feo.55Mno.45Ti03 under a dc bias field H at temperature 
T. We first specify, in a manner explained in detail be- 
low, sets of (T, H) , at which the characteristic relaxation 
time of the system, t{T;H), coincides with the inverse 
of the ac field frequency uj, i.e., 1/uj = t{T; H). We then 
examine whether or not these t(T; H) obey a dynamical 
critical behavior represented by 



t{T-H) ^ hiiJL^Y - to\T/TgiH) - ir 



(1) 



at T > Tg{H). Here we suppose that the SG replica sym- 
metry breaking (RSB) phase predicted to appear below 
the AT line, Tg{H), by the mean-field theory is accom- 
panied with a high-temperature disordered phase which 
exhibits critical divergences of the correlation length 
^c{T;H) - Lo\T/Tg{H) - ll'" and of the correlation 
time given by the above equation, where to and Lq are 
respectively microscopic units of time and length, and z 
and v are critical exponents. It turns out that the data 
t(T; H) are not compatible with the expected dynamical 
critical behavior except for the case with H — 0. 

In the droplet theory, on the other hand, the so-called 
magnetic correlation length, ^h, is introduced. It spec- 



ifies the behavior of droplet excitations with a linear 
size L in a SG state under a field H Q. For droplets 
with L < their behavior is governed by the SG 
stiffness free energy T{T){L/ Lq)^ , while for those with 
L > S.H their behavior is dominated by the Zeeman en- 
ergy ^/q^HiL/Laf/^, where T(T) is the SG stiffness 
modulus, 9 the stiffness exponent, gEA(T) the EA order 
parameter, and d the spatial dimension. Explicitly, at T 
less than Tc[= Tg{Q)\ which is a unique critical point of 
the system, is written as 



T(T) 



i?V9^I(T) 



(1 - r/Tc)"-"' 

H 



(2) 

where 5 = 2/{d — 29). In the last expression the 
temperature dependence of T / ^gsA is represented by 
(1 - T/Tc)"'^". AtT < Tc, acff = 9iy - (3/2 is expected, 
where (3 is the critical exponent of qEA- In the present 
analysis on the ac susceptibility measurements of fre- 
quency LO under H, we identify S^h with LT{t ~ 1/^^)1 
which is the mean size of droplets that can respond to 
the ac field of frequency uj at temperature T, i.e.. 



(3) 



Furthermore it is considered that Lxit) has the same 
functional form as that of the growth law Rrit) of the 
SG correlation length which grows with time t after the 
system is kept at a constant temperature T under H = 0. 
Explicitly, we adopt here an algebraic growth law 



Lrit) ^ Rrit) ^ Loit/to)'^'^^ 



(4) 



which is commonly observed in numerical simulations [l^ 
0, ^3 as well as in experiments 0, . Interestingly, 
the present data r(T; H) for relatively large H and so 
relatively low T turn out to be consistent with the sce- 
nario represented by Eqs. JSJ, and iQ. We interpret 
this result as evidence for the presence of the magnetic 
correlation length S^h predicted by the droplet theory. 



2 




10 20 30 40 



FIG. 1: (Color online) ac susceptibility vs temperature for 
bias fields H ^ 0, 100, 300, 1000, 3000, 6000, 10 000, 15 000, 
and 20 000 Oe. 



M (measured with a slow cooling rate) coincide with 
each other, as Xc,(T,i?) = '-^^^ « ^Mili = 

order to probe the linear response. 
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FIG. 2: (color online) (a) MzFc/h, h'hc/h and x'{T,H;uj) vs 
temperature in zero bias field, (b) Xcq and x' vs temperature 
in a bias field of H = 1000 Oe. Dots indicate the condition of 
Eq. @. 



We measure the ac susceptibility, xi'^iT.H) = 
x'(w,T,i7) + ix"{uj,T,H), of the Ising spin glass 
Feo.ssMno. 45X103 in the frequency range of 0.001 s < 
1/u! <1 s. All measurements are performed on a MPMS- 
XL squid magnetometer equipped with the low-field op- 
tion. The susceptibilities x' and x" measured under dif- 
ferent bias fields are shown in Fig.^ while shown in Fig. [21 
are the dc and ac susceptibilities of different frequencies 
in zero bias field and in a bias field oi H = 1000 Oe. 
The equilibrium susceptibility Xeq(T,H) can be deter- 
mined in the temperature range within which the dc field- 
cooled (FC) and zero-field-cooled (ZFC) magnetization 



FIG. 3: (Color online) [xcq — x'(^)]/Xcq time on a dou- 
ble logarithmic scale. The lines are fits to a power law for 
H —0 and 300 Oe or to a stretched-exponential for H =1000, 
2000 and 3000 Oe. Inset: the sets of {T,H) determined by 
Eq. (gj forr(r;iy) equal to 9.4- 10"\ 3.1-10"\ 9.4■10"^ 3.1- 
10~^ 9.4 ■ 10"^ 2.8 ■ 10"^ 9.4 ■ lO^* s (left to right). 



As indicated by circles in Fig. 13 we identify the char- 
acteristic relaxation time t{T;H) (= l/ui) by the condi- 
tion [l3, 



x\uj,T,H) = 0.98xcq{T,H 



(5) 



To explain our idea behind this condition we present in 
Fig-ElXoq ~ x'(i^)i which is proportional to the dynamic 
spin correlation function q{t) with t = 1/lo, for different 
bias fields at T = 20 K. It has been found to follow the 
empirical formula q(t) = Ci^^e^^*/'^*^-^)''' at T above 
Tc and a pure power law below Tc for H = both in 
numerical simulations on the EA Ising model0] and in 
experiments on Fen »^Mnn »^TiO^[l9|. The results are in- 
terpreted to indicate that the distribution of relaxation 
times is bounded from above by a certain value around 
T*(T) at T > Tc, while r*(T) is infinite at T < Tc. 
Although our experimental timewindow for a fixed set of 
(T, H) is rather limited, Fig.Oindicates that it is possible 
to bring the Ising SG to equilibrium by applying a mag- 
netic field. In particular, for the set of (T, H) for which 
Eq. lO is found to be satisfied, the corresponding q{t) 
exhibits a stretched-exponential form and so the corre- 
sponding T*{T; H) is definitely finite. We have therefore 
simply introduced the condition of Eq. Q to specify the 
upper bound of relaxation times without fitting our q(t) 
to a stretched-exponential form explicitly. We consider 
that the present method for specifying t{T] H) is system- 
atic and preferable to other methods used previously ji^. 
It is also noted that the condition is satisfied in the tem- 
perature range where we can determine Xcq as seen in 
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Fig. [3 implying that t(T; H) thus specified in fact asso- 
ciates with certain relaxation processes in equilibrium. 
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Tg(H) = T/[1 +(T/y(-""">] 

FIG. 4: (Color online) Test for critical dynamics assum- 
ing that the spin glass transition temperature T^{H) — 
T{t;H)/[1 + {r/to)'^/'""] [Eq. ||TJ] changes with r for each 
H as explained in the text. The inset shows r vs T(r; H = 
0)/Tc - 1 with Tc = Tg{H = 0) = 22.3 K on a log-log scale. 

In the inset of Fig. |2| we show the H — T rela- 
tions which yield common values of t{T;H). They are 
roughly consistent with the AT-lines, i.e., of the form 
H [1 - T(r;7J)/r(r;0)]3/^ However, such appar- 
ent AT-lines dependent on measuring time scales r are 
by no means a proof of an equilibrium SG phase un- 
der a field We require the dynamical critical slow- 
ing down represented by Eq. ^ to hold at T > Tg{H) 
where Tg(H) is an assumed critical temperature under 
a field H. As shown in the inset of Fig. ^ t(T; 0) in 
our timewindow are in fact fitted to this expression with 
Tc = Tg(0) w 22.3 K, io « 3 • 10"" s, and zv w 11. The 
value of zv is in accordance with the previous values ob- 
tained for the Ising spin glass Feo.soMno.soTiOa 0- This 
result gives evidence for an equilibrium SG phase tran- 
sition in zero field. The main frame of Fig. 0] is, on the 
other hand, the result of an attempt to see if the system 
also exhibits critical slowing down in small dc fields as- 
suming that to remains the same as in zero field. The 
best fit to Eq. under this constraint is found with 
ziy w 22 for all H > 0. To demonstrate the quality of 
the fit we show a. H — T diagram whose T-axis is Tg{H) 
calculated for each t(T; H) using the obtained value of 
zv. A unique Tg{H) is found only for H — 100 Oe (and 
H = Q) and the data in iJ > 100 Oe exhibit systematic 
dispersion, giving evidence against critical dynamics for 
H > 100 Oe. Even for H = 100 Oe, however, due to 
the large value of ziy « 22 we consider that the relatively 
good fit to Eq. has no physical meaning. We also note 
that if both to and zv are adjusted in the fitting, the re- 
sultant to takes unphysically large values. We therefore 
conclude that, under a magnetic field, the system does 
not exhibit a phase transition which is accompanied with 
critical slowing down. 

Next, let us examine the data t{T; H) for relatively 




a^„ln(1-T/T^)-ln(H) 

FIG. 5: (Color online) The scaling plots of {T/Tc)\nT{T- H) 
vs Oeffln(l - T/Tc) - \nH with a^s = 0.25. Open (filled) 
symbols mark temperatures > (<) O.TTc. 



large H based on the droplet picture, namely, by regard- 
ing r(T; H) as a relaxation time of the largest SG clusters 
of a mean size £,h{T; H), which is determined by Eq. Q 
combined with Eqs. Q and Q). For this purpose we 
plot (r/Tc)ln(T/to) for H > 5000 Oe as a function of 
aoffln(l — T/Tc) — InH, thereby adjusting aoff but keep- 
ing to fixed to the value obtained above. As shown in 
Fig. the best collapse of the data at T/Tc < 0.7 is 
obtained with aoff « 0.25. The slope of the fit, which is 
equal to b/S, gives b w 0.11, where 9 = 0.2 (and 
so 6 —0.77) are used. Here we examine T{t;H) in a 
rather narrow temperature range (the lower bound of 
our observation of r is 0.55Tc), since the approximation 
T~/%/9ea ~ (1 ~ T/TcY"" used to derive the last expres- 
sion of Eq. (121) with a constant Ccff is not expected to work 
in a wider temperature range. In fact when we analyze 
the r(T; i?) data obtained for Feo.sMuo.sTiOa by Matts- 
son et al. |3 by the present method, we obtain acff ~ 
and b k. 0.09 for their data at T < 0.5Tc, and aoff w 0.2 
and b « 0.09 for their data at 0.5Tc < T < 0.7Tc. 
The circumstances become more subtle at T closer to 
Tc, since the critical behavior of T(T) and qEA{T) has 
not been well established yet. We can only mention that 
the same analysis on our data t{T; H) at T < 0.95rc 
and for H < 5000 Oe yields aoff ~ 0.45 and b « 0.13. 
This strongly implies that the critical exponent, Oefi at 
T = To, is positive. In spite of such subtlety concerning 
with the temperature dependences of T(r) and q^AiT), 
the results shown in Fig. |S1 are sufficient for us to con- 
clude that the system is in the paramagnetic phase, with 
and r(T; H) being the upper bounds for the SG cor- 
relation length and relaxation time, respectively. 

The values of the exponents extracted above from 
t{T; H) at 0.55Tc < T < 0.7ro can be compared with 
those obtained in the simulation on the Ising EA model 
in the corresponding temperature range; b « 0.16 [l3 | 
and aoff ~ 0.14, where the latter value is extracted from 
S^H which is obtained and denoted as ItH~^ in Q. These 
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figures are in turn compatible with tiie ZFC magneti- 
zation of Fen rMuo RTiO?; observed in heating processes 
with intermittent stops [l5j as well as with the crossover 
scenario for the AT-like transition observed also in 
Feo.sMno.sTiOa 0]. Here we emphasize that the values 
of b extracted from the simulation and the experiments, 
over more than 10 decades difference in time scales in 
units of to, agree with each other reasonably well. 

One more comment is in order on the growth law of the 
SG correlation at a constant T. In the original droplet 
theory instead of Eq. (0J, a logarithmic form of 



ported by the NAREGI Nanoscience Project, from the 
Ministry of Education, Culture, Sports, Science, and 
Technology. 



Lrit) - Lo 
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1/-4, 
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has been proposed, where A is the characteristic energy 
scale of the frcc-encrgy barrier. For Fco.soMno.soTiOs 
Eq. ((nj has been applied to its various phenomena 
with the resultant exponent ip ranging from 0.03 to 
1-9 0,111 When Eq. © is replaced by Eq. ® 

with a T-independent A in the present analysis, we ob- 
tain -0 « 0.6 and the same as obtained above from 
the t{T;H) data at O.SSTc < T < O.TTc. We consider, 
however, that the power-law growth is superior to the 
logarithmic growth in the sense that the values of b ex- 
tracted from various phenomena arc less dispersive than 
those of ip . The power-law growth implies that the 
free-energy barrier of droplets excitations of a size L is 
proportional to InL which becomes smaller than , the 
free-energy gap of the corresponding droplet excitations, 
for a sufficiently large L, say L*. Therefore, our conclu- 
sion that the power-law is a more plausible description 
of various SG glassy dynamics mentioned above implies 
that the SG short-range of less than L* is involved in 
such slow-dynamical phenomena observed even by lab- 
oratory experiments. In particular, Xcq and x'{^) ana- 
lyzed in the present work involve length scales S,c{T) and 
Lt{1/uj) or which arc much shorter than L* so that 
they can be regarded as SG equilibrium properties (see 
[lll | for non-equilibrium phenomena associated with sim- 
ilar length scales much less than L*). We believe that, in 
the field of glassy dynamics, proper understanding of the 
length and time scales of phenomena we observe is very 
important, though it is often not an easy task. 

To summarize, we have found experimental evidence 
against the existence of an equilibrium phase transition 
in Ising spin glasses under a bias magnetic field, i.e. evi- 
dence against the AT phase transition. From the present 
analysis we learn that one has also to be careful to 
draw conclusions about the equilibrium phase diagram 
of Heisenberg spin glass under a field; the occurrence of 
the irreversibility [ij alone may not be evidence for the 
presence of an equilibrium phase transition. 
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